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ABSTRACT 

By a systematic method we construct general classes of exact and selfconsistent ax- 
isymmetric MHD solutions describing flows which originate at the near environment 
of a central gravitating astrophysical object. The unifying scheme contains two large 
groups of exact MHD outflow models, (I) meridionally self-similar ones with spherical 
critical surfaces and (II) radially self-similar models with conical critical surfaces. The 
classification includes known polytropic models, such as the classical Parker model 
of a stellar wind and the Blandford and Payne (1982) model of a disk-wind; it also 
contains nonpolytropic models, such as those of winds/jets in Sauty and Tsinganos 
(1994), Lima et al (1996) and Trussoni et al (1997). Besides the unification of these 
known cases under a common scheme, several new classes emerge and some are briefly 
analysed; they could be explored for a further understanding of the physical properties 
of MHD outflows from various magnetized and rotating astrophysical objects in stellar 
or galactic systems. 

Key words: MHD - plasmas - solar wind - stars: mass loss, atmosphere - ISM: jets 
and outflows - galaxies: jets 



1 INTRODUCTION 

A widespread phenomenon is astrophysics is the outflow 
of plasma from the environment of stellar or galactic ob- 
jects, either in the form of a noncollimated wind (Parker 
1958, Feldman et al 1996), or, in the form of collimated jets 
(Blandford & Rees 1974, Biretta 1996, Ferrari et al, 1996). 
These outflows not only occur around typical stars and the 
nuclei of many radio galaxies and quasars, but they are also 
associated with young stars, older mass losing stars and 
planetary nebulae, symbiotic stars, black hole X-ray tran- 
sients, low- and high-mass X-ray binaries and cataclysmic 
variables (for recent reviews see e.g., Ray 1996, Kafatos 
1996, Mirabel & Rodriguez 1996, Livio 1997). Even for the 
two spectacular rings seen with the HST in SN1987A, it has 
been proposed that they may be inscribed by two precessing 
jets from an object similar to SS433 on a hourglass-shaped 
cavity which is created by nonuniform winds of the pro- 
genitor star (Burderi and King, 1995, Burrows et al 1995). 
Also recently, in the well known long jet of the distant radio 
galaxy NGC 6251 an about 10 3 light-year-wide warped dust 
disk perpendicular to the main jet's axis has been observed 
by HST to surround and reflect UV light from the bright 
core of the galaxy which probably hosts a black hole (Crane 
& Vernet 1997). 



Nevertheless, despite their abundance the questions of 
the formation, acceleration and propagation of nonuniform 
winds and jets have not been fully resolved. One of the main 
difficulties in dealing with the theoretical problem posed by 
cosmical outflows is that their dynamics needs to be de- 
scribed - even to lowest order - by the highly intractable 
set of the MHD equations. As is well known, this is a non- 
linear system of partial differential equations with several 
critical points, etc, and only very few classes of solutions 
are available for axisymmetric systems obtained by assum- 
ing a separation of variables in several key functions. This 
hypothesis allows an analysis in a 2-D geometry of the full 
MHD equations which reduce then to a system of ordinary 
differential equations. The basis of the self-similarity treat- 
ment is the prescription of a scaling law in the variables as a 
function of one of the coordinates. The choice of the scaling 
variable depends on the specific astrophysical problem. 

In spherical coordinates (r,6,<f>), a first broad class 
for describing outflows are the so-called meridionally self- 
similar MHD models. Parker's (1958) classical modeling of 
the spherically symmetric polytropic solar wind is the sim- 
plest member of this class. A new class of such type of mod- 
els for describing magnetized and rotating MHD outflows 
from a central gravitating object has also been examined 
(Sauty & Tsinganos 1994 (henceforth ST94), Lima et al, 
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1996, Trussoni et al 1997). For example, an energetic crite- 
rion for the transition of an asymptotically conical outflow 
from an inefficient magnetic rotator to an asymptotically 
cylindrical outflow from an efficient magnetic rotator was 
derived. In the present paper, it will be shown that this 
special class of meridionally self-similar solutions is one of 
the simplest possible meridionally self-similar models. Fur- 
thermore, a new interesting member of this class of radially 
self-similar MHD models will be briefly sketched. 

A second broad class of solutions contains the radially 
self-similar MHD models. Bardeen & Berger (1978) pre- 
sented the first such models in the context of hydrodynamic 
and polytropic galactic winds. Nevertheless, their generali- 
sation to a cold magnetized plasma by Blandford & Payne 
1982 (henceforth BP82), remains widely known because of 
their success in showing for the first time that astrophysical 
jets can be accelerated magnetocentrifugally from a Keple- 
rian accretion disk, if the poloidal fieldlines are inclined by 
an angle of 60°, or less, to the disk midplane (but see also, 
Cao 1997). A further extension has been presented by Con- 
topoulos & Lovelace (1994) for a hot plasma with a more 
general parametrization of the magnetic flux on the disc, 
while these models form the basis of several investigations 
of accretion-ejection flows from stars and AGN (Konigl 1989; 
Ferreira & Pelletier 1995; Ferreira 1997; Li 1995, 1996). In 
this paper it will be shown that this special class of radi- 
ally self-similar solutions is one of the simplest possible such 
models. Furthermore, a new interesting member of the ra- 
dially self-similar MHD models will be sketched. 

In subsection 2.1 we use a simple theorem in order 
to construct several classes of meridionally selfsimilar so- 
lutions and the resulting cases are then summarized in Ta- 
bles 1 and 2. The general method is next applied in sub- 
section 2.2 to a step by step construction of a new model 
for collimated outflows which is also briefly sketched there. 
In section 3 the other remaining possibility in spherical co- 
ordinates, i.e., radial self similarity is taken up. The result- 
ing cases are summarized in Table 3 while a new model is 
also briefly sketched which gives asymptotically cylindrical, 
paraboloidal and conical streamlines. Finally, the results are 
summarized in Sec. 4. 



2 MERIDIONALLY SELFSIMILAR 
MHD OUTFLOWS 

Consider the steady (d/dt = 0) hydromagnetic equations. 
They consist of a set of eight coupled, nonlinear, partial 
differential equations expressing momentum, magnetic and 
mass flux conservation, together with Faraday's law of in- 
duction in the ideal MHD limit, 



(V x B) x B 



47T 



VP - pVV , 



V-B = 0, V-(pV)=0, Vx(VxB)=0. 



(1) 



(2) 



B, V, — VV = —V (— QM/r) denote the magnetic, velocity 
and external gravity fields, respectively, while p and P the 
gas density and pressure. With axisymmetry {d/d<f> = 0), we 
may introduce the magnetic flux function A, such that three 
free integrals exist for the total specific angular momentum 



carried by the flow and the magnetic field, L(A), the coro- 
tation angular velocity of each streamline at the base of the 
flow, Q(A) and the ratio of the mass and magnetic fluxes, 
^a(A) (Tsinganos 1982). In terms of these integrals and 
the square of the poloidal Alvfen Mach number (or Alfven 
number), 



47rp ' 



(3) 



the magnetic field and bulk flow speed are given in spherical 
coordinates (r, 0, <f) 



B = Vx 



A(r,6)<f> L^a - r 2 sin 2 0fi*, 



r sin 0(1 - M 2 ) 



4-zrp r sin ( 



+ 



r_ sin 2 OQ, - LM 
r sin 0(1 - M 2 ) 



(4) 



(5) 



To construct classes of exact solutions, we shall make two 
crucial assumptions: 

(i) that the Alfven number M is some function of the 
dimensionless radial distance R — r/r+, 



M = M(R) , 
and 



(6) 



(ii) that the poloidal velocity and magnetic fields have a 
dipolar angular dependence, 



A = 



r IB, 



A (a) , 



R 2 



G 2 (R) 



(7) 



By choosing G (R — 1) = 1 at the Alfven transition R = 1, 
G(R) evidently measures the cylindrical distance vo to the 
polar axis of each fieldline labeled by a, normalized to its 
cylindrical distance vj a at the Alfven point, G (R) = vj /vj a . 
For a smooth crossing of the Alfven sphere R = 1 [r = 
r+,6 = 9 a (a)], the free integrals L and fl are related by 

^ = ™!{A) = r 2 sin 2 a (a) = rla . (8) 

Therefore, the second assumption is equivalent with the 
statement that at the Alfven surface the cylindrical distance 
zu a of each magnetic flux surface a — const is simply pro- 
portional to \[a.. 

Note also that the gravitational potential can be expressed 
in terms of the escape speed V eS c at the Alfven radius r*, 



V 



2R 



Vesc 



2QM 



Instead of using the three free functions of a, (.4 , ^ a , 
we found it more convenient to work instead with the three 
dimensionless functions of a, (g\ , <?2 , 53), 



9i (a) 



A 2 da , 



32 (a) = / 2 *ida 



93 (a) = 



47T/9* 



(9) 
(10) 
(11) 
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Substituting for II in Eq. (|l^) it follows, 
- fcg-i - faag 2 + f 5 g 2 - fiag x - f 6 g 1 + f 4 gi + f = ,(23) 



Also, we shall indicate by II the total pressure in units of the 
magnetic pressure at the Alfven surface on the polar axis, 
BJ/8tt = p*K 2 /2, 



8tt 



such that, 



B- 



P = 



8tt 



n + /iffi + fcag-t + f 3 ag 2 



(12) 



The functions fi(R), i ~ 1,2,3 are given in Appendix A 
while all starred quantities refer to their respective values 
at the polar Alfven point (R — 1 ,a — 0). Hence, 



A (a = 0) = 1 , *a (a = 0) = ^471-/9* 
or, 

g[ (a = 0) = 1 , g 3 (a = 0) = 1 . 



(13) 



With assumptions (i)-(ii) and in this notation, the f- 
and 9— components of the momentum equation become, 

dU (R,9) , ( F \ > 

F \ 



(14) 



2cot0 ( faagx + f 5 ag 2 



(15) 



dll (R, 9) 
80 

Next, by using a instead of 9 as an independent variable, we 
may transform from pair of the independent variables (R , 6) 
to pair of the independent variables (R,a). With the follow- 
ing elementary relations valid for any differentiable function 



d$(fl,6>) 

M 

d$(R,9) 



d$(R,a) Fd$(R,a) 
+ a- 



— 2a cot 



OR R da 

&f>(R,a) 



da 



(16) 



(17) 



we may transform Eqs. (14), (|15j) into the following two 
equations: 

dU (a, R) 



h9i + h a 9i + hag 2 + f 9 g 3 ■ 



(18) 



(19) 



(20) 



da 

dU (a, R) 
dR 

By integrating Eq. ( jli| ) we get II = f4g± + /5O2 + fo where 
fo is an arbitrary function of R. From Eq. ( |l2| ) the pressure 
is 

87, 

or, 

8tt 

where P and Y are the (1 x 7) matrices, 

p = [/o h h h h h o] , 

and 

Y = [ Yi Y 2 Y 3 Yi Y 5 Y 6 Y 7 ] = 



3^ f tit 
j~~ ( fagi + h92 + fo + figi + fcagx + f 3 ag 2 

-Ypt , 



(21) 



1 9i 9i agi 92 ag 2 g 3 



(22) 



an expression of the form 

X 7 {R) Y 7 (a) + X 6 (R) Y 6 (a)+--+Xi (R) Yi (a) = 
or , Yxt = 

with X the (1x7) matrix 
X = [ Xi X 2 X 3 X4 X5 Xg X7 ] = 
fo fi - / 6 - / 7 fs —fa - h 



(24) 



(25) 



2.1 Systematic construction of classes of 

meridionally selfsimilar MHD outflows 

In (Vlahakis & Tsinganos 1997) the following simple theo- 
rem was proved: 

Theorem: If F n (a), Yi(a), Xi(R), i — 1,2, ■•■,n are arbi- 
trary functions of the independent variables a and R and 

F n (a) = Fj (a) Xi (R) H +Y n (a) X n (R) , (26) 

then, there exist constants Ci, c 2 , . . . , c n such that, 

F n (a) = ciYj (a) + c 2 Y 2 (a) + ■ ■ ■ + c n Y n (a) . (27) 

Consider then a relation of the form, 

X n (R) Y n (a) + ■ ■ ■ + Xi (R) Yi (a) = . (28) 

Regarding the first term of the sum there are evidently only 
two possibilities. Either, 

(i) X n (R) — for every R, in which case (indicated by 
the digit "0") we have 

I„-i (R) Y„_! (a) + • • • + Xi (R) Y a (a) = , 

or, 

(ii) Xn (R) ^ 0, in which case (indicated by the digit "1") 
we have 



Y n {a) 



X! (R) 



Yi(a) 



-Yn-i (a) 



X n (R) * 1 X n (R) 

Then, according to the theorem stated in the be- 
ginning of this section, there are constants = 
1, 2, • • • , n - 1 such that Y n (a) = m!™' Yi (a). This 

gives a condition between the functions of a. Substitut- 
ing this in the initial sum we find: 



+ 



A„_! (R) + It^Xn (R) 

A„_ 2 (R) + ^ n %X n (R) 
Xi (R) + ^X n (R) 



Y„_! (a) + 
Y„_ 2 (a) + ■ 
Yi (a) = . 



(29) 



Hence, in both cases we find a sum with n — 1 terms. Fol- 
lowing this algorithm at the end we '11 have only one term. 
Since for each product we have the above two possibilities, 
totally we obtain 2" cases. Each of them corresponds to a 
set "xx- • -xx" with x — 1,0 (n digits). The number of "1" 
digits is the number of conditions between functions of a 
while the number of "0" digits is the number of conditions 
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Table 1. Meridionally Selfsimilar Models 

Case gi (a) 92(0) 9'i( a ) constants 

(1) a X 2 a 1 + Sa 

(2) a ^a + iia^/e 1 + 8a + a8 a e e^0,l,u^0 

(3) a ^a + ualna 1 + 8a + fiSga In o ^7^0 

(4) aoe°o Ae=o 1 + <5ae "0 + — 1^ 

(5) Sfl. I -2 1 (-2 l) f I -2 1 |e 1 + 5 I _2_ _ ! 16 + I _2L _ 1 |e-l 6 ^0,1 

(6) -a ln|-^-l| ^ l n I ^ 1 | 1 + S In | -2L - 1 | + u—^ . 

(7) /t ln^+Ca 1 + 5( Q ,_ Q , ) + /Lt ,5 l n -^ a ^ 

(8) eTi^T (=j)' Aia' + Aaa- 1 1 + Si (a« - ag) + <5 2 (a- 1 - ag" 1 ) ^0,1 

(9) t^I- In Ai In + ^ 1 + <5 a In SL + 82 ( - - — ) 



between functions of ii. 

Now, following this method from Eq. we get 2 7 

solutions. Each of them corresponds to a set "xxxxxxx" with 
x either 1, or, 0. Of those numbers: 

(i) The first digit is always "1" (because X7 7^ 0). 

(ii) The last digit is always "0" (because Y\ 7^ 0). 

(iii) Since A 1 7^ 0, it follows that g[ 7^ and thus 31 
cannot be a constant. Hence, the function Y% = gi cannot 
be proportional to Y\ and therefore all numbers always have 
"00" at the end. 

(iv) We have totally six unknown functions: the three 
functions of R, (G,M,fo) and the three functions of a, 
(31,32, 33)- On the other hand, the number of conditions 
between the functions of R (their number is equal to the 
number of digits "0") and the functions of a (their number 
equals to the number of digits "1") in each one of the sets 
"xxxxxxx" is seven. It follows that the system of (G, M, fo) 
and (gi, 32, 33) is overdetermined. Note however that since 
the forms of the functions Xi(R) are much more compli- 
cated than the forms of the functions Yi (a), we choose sets 
"xxxxxxx" with at most three "0's" because in the case of 
4 or more "0's" we have correspondingly 4 or more relations 
between the 3 functions of 7?, which in general overdeter- 
mines the system of (G, M, fo). We then shift the problem 
of overdetermination of the problem to the set of the 3 func- 
tions of a, (<7i, g2, 33) which need to satisfy 4 relations. In 
this system however, it is possible to choose the constants 
/ij such that a consistent solution for the functions of a 
can be finally constructed. 

Altogether, then and with these considerations in mind, 
from the 2 7 = 128 possible cases we end up with only five: 
1011100, 1101100, 1110100, 1111000, 1111100. For each of 
one of those sets we can solve the system for 31 , 32 , 33 , as it 
is shown in the example of the next section. 
From a different perspective, 31(a), 32(a), 33(a) are vectors 
in a 3D a-space with basis vectors [«i(a) , fi2(a) , 113(a)]. 
This space contains all vectors 3i(a), i = 1,2,3 subject 
to the #-selfsimilarity constraint manifested by Eq. (E3[), 
i.e., that for a given such set gi(a), i = 1,2,3, the vectors 
1,03! (a) i a 32(a) and 31(a) also belong to the same space. 
Each of the resulting functions gi(a), i = 1,2,3 are then a 
linear combination of the basis vectors Ui(a), W2(a), 113(0). 
In the following, we choose ui — 1, U2 — 31(a). All such 
sets of basis vectors give all possible meridionally selfsimi- 
lar solutions. Therefore, collecting all possibilities, we end 



up with the classes of solutions shown in Table 1. Note 
that in the last three cases A (a = 0) 7^ 1, but one can 
say that the starred quantities refer to values at the point 
R = 1, a = ao < 1. 

In all nine cases of Table 1, from Eqs. (§), @, §H§ 
we may find easily the forms of the free integrals from the 
relations, 



.4 



33 ■ 



L = r i ,V i ,a\ 



(30) 



(31) 



while by substituting 31,32,33 in Eqs. (J20|), (|23|), the cor- 
responding ordinary differential equations for the jet radius 
G(R), Alfven number M(R) and pressure component fo(R) 
are found from the i?-relations, as it is illustrated in the fol- 
lowing section. 

From the perspective of the a-space, in each one of the cases 
of Table 1 there exists a 3 x 7 matrix K such that 



(32) 



Y = [ Ml U2 U3 ] K , 

so that from Eq. (24), 

[ Ml U2 U 3 } KX) = . 

If Ui are linearly independent then 
KX* = . 

These three equations are the ordinary differential equations 
for the functions of R in each model while the pressure is, 



P = ft [ ui U2 u s ] Kpt = ft (p . 
Sit Sit 

where 



P131 + P2M3) , 



KP 



t 



[Po Pi P2 



it 



The first two cases of Table 1 are of some interest. The first, 
is a degenerate one with U3 = and the following form of 
the free integrals: 



A = ^ 



a, <f A = \747rp* (1 + 5a) , 



(33) 



This is a special case of the more general following case (2) 
for u — (and £ = A 2 ) and has been studied in detail in 
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Table 2. Mcridionally Sclfsimilar Radial Models 



Case 


91(a) 


92(a) 


93(a) 


(1) 


- In | 1 - a | 





1 


(2) 


^ / i^<to-ln| 1-a | 


A 2 — 

e 


1 + <5a e 


(3) 


Mi In | l-a\ + l M2$ ^da 


A In a 


5i + <52 In a 


(4) 


gi(a) ^ fii In | 1 - a | + 1x2 





Sg[ (1 — a) + A 


(5) 


p, In | 1 — a | 


92(a) ^ (p, pi hi a + p2,pia^ + ^3) 


5 



ST94 and Trussoni et al (1997). It is the single case where we 
have only two conditions between the functions of 7?, so that 
the third relation between the unknown functions G, M, fo 
is freely chosen. In Trussoni et al (1997) this corresponds to 
an a priori specification of the shape of the poloidal stream- 
lines, while in ST94 in an a priori imposed relationship be- 
tween the spherically and nonspherically symmetric compo- 
nents of the pressure. This last case leads to a generalized 
polytropic-type relation between pressure and density of the 
form, 



P(a,R) p(a,R) 
= function of 



(34) 



P(0,R) P (0,R) 

As a result, a Bernoulli- type constant exists and, among 
others, this constant gives a quantitative criterion for the 
transition of an asymptotically conical wind from an ineffi- 
cient magnetic rotator to an asymptotically cylindrical jet 
from an efficient magnetic rotator. 

The second case with e / 0,1,^ ^ has U2 = a, U3 — a € . 
The corresponding form of the free integrals is : 



A = ^-a, *a = v / 4 : irp+(l+Sa + fj,S at) , 



Q 



(35) 



This is a new case which emerged from the present system- 
atic construction. The corresponding differential equations 
are derived in detail in the example of the next section where 
the solution is briefly analysed. 

In the special configuration with G = R a = sin 2 9, the 
field and stream lines on the poloidal plane are radial and 
we find five cases shown in Table 2. 

The first case is a degenerate one, wherein there is only 
one condition between the unknown functions M(R) , fo(R)- 
Thus, a second relation between M(R) — fo(R) can be 
imposed a priori, for example, a polytropic relation be- 
tween pressure and density. This last possibility leads pre- 
cisely to Parker's (1963) classical solar wind solution with 
a radial and nonrotating outflow. All other cases (2)- (5) 
are non-degenerate, i.e., there are two relations between 
M(R) - f (R). 

The second case has been analysed in detail in Lima et al 
(1996) and corresponds to a radial but heliolatitudinally de- 
pended outflow. In addition fi = —1, e = 1 this case co- 
incides with (1) in Table 1 for radial poloidal streamlines. 
Note that a common feature of all rotating cases with radial 
stream lines on the poloidal plane is that they cannot be ex- 
tended in all the poloidal plane, for sufficiently fast magnetic 
rotators. For example, in the model of Lima et al. (1996) the 
pressure becomes negative at some colatitude 9 max , for large 
values of rotation. This is basically due to the fact that with 
the poloidal magnetic field dropping like 1/J? 2 and the az- 



imuthal field dropping like 1 / R, the magnetic pressure drops 
like 1/R 4 and by itself alone cannot balance the magnetic 
tension which drops like 1/R 3 ; a strong pressure gradient is 
then needed from the pole towards the equator to balance 
the magnetic pinching. In fast magnetic rotators this pres- 
sure gradient is so strong that it leads to negative values of 
the pressure at angles 6 > 9 ma x- A collimated outflow with 
uniform asymptotic conditions is the only way left for an 
everywhere valid outflow from an efficient magnetic rotator 
(Heyvaerts & Norman 1989, ST94). 



2.2 Example of a new model for a meridionally 
self-similar MHD outflow 

Let us illustrate the previous construction with the example 
1101100 obtained from the present case with n = 7. This 
number means the following: 

Since the first digit is 1, there are six constants p\ ,i = 
1, 2, • • • , 6 such that the following relation holds between the 
functions Yi(a), i=l,2, ..7, 



Y 7 = ^2 Y i , (a-relation-1) . 



(36) 



Substituting this expression of Y7 in the initial relation Eq. 
(24) between the functions (Xi, Yf), i=l,..7, we obtain 



(x + H ( pX 7 ^ Y 6 + (x 5 + n ( pX 7 ^ Yb + • • • + 
(Xi +a4 7) X 7 ) Yi =0. 



(37) 



Now the second digit is again 1 and thus there are five con- 



stants p!f\i = 1, 2, • • • , 5 such that 

5 

Ya = Yl Yl ' (a-relation-2) , 

i=l 

while substituting this relation in Eq. (37) we obtain, 
X 5 + p^Xr^j + M < 6) (x 6 + p { PX 

(Xi + 4Px 7 ^j + PP (x + PPX 7 J 
The third digit is and hence 

(x s + PPxPj + ^ (x 6 + PpxPj = 

(ii-relation-1) 



(38) 



Y 5 + ■ ■ • 
Yi =0. 



(39) 



(40) 
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a relation between the functions of R. With the help of Eq. 
(40), Eq. (39) now reduces to, 



Y t = 0. (41) 



The fourth digit is 1 and thus there are three constants 
nP ,i = l,2,3 such that 



Ya = /4 4 ' Yi , (a- relation- 3) . 
Substituting this relation in Eq. (41) we obtain 



(42) 



+ 



,( 4 ) 



(Xa + p^Xr) + fJ>f ] (x e + Ms 7 '^) } Yi = 



(43) 



The fifth digit is 1 and there are two constants /4 , i = 1, 2 
such that 



y 3 = /4 3) yi + nt'Yz , (a-relation-4) . 



,0) 



(44) 



Substituting this in Eq. (43) we find a relation involving Y\ 
and Y-i . Finally, we must put equal to zero the multipliers of 
Yi,Yi in this relation because the two remaining digits are 
0. So we have 



(Xi + uPXr) 

(x 4 + ^ ] X 7 ^ + M f (x e + $ ] x 7 \ 
X 3 +^X r )+^ (Xe + uPXr) 



,(3) 



(4) 



X 4 + fi { Px 7 ) + fi> ( X 6 + n K 6 "X 7 



,(7) 



+ 



= 0, 



(45) 



for i = 1, 2(P-relations-2,3) . 

These last two equations together with Eq. (40) are the three 
equations for the functions of R. On the other hand, Eq. 
(36), Eq. (38), Eq. (42) and Eq. (44) are four relations among 
the three functions of a. These relations of the functions of 
q [Eqs. (44), (42), (38), (36)] are equivalent to the system: 



Y 3 =c 1 Y 1 + c 2 Y 2 

Y4, = C3F1 + C4Y2 

Y a = c 5 Yi + c 6 y 2 + C7F5 

Yr = c 8 Yi + C9Y2 + C10Y5 



g x = ci + c 2 gi 
Qffi = C3 + c 4 gi 

«ff2 = C5 + C 6 ffl + C 7 #2 

ff3 = c 8 + c 9 gi + ciog2 



Note that we renamed the constants and also used Eq. (22). 
From the first, if C2 7^ it follows that gi = — C1/C2 + ce C2< *. 
Then, from the second c = and hence gi — —c\jci. But 
gi cannot be a constant. Thus, C2 = while the first two 
equations combined with Eq. ( p^ ) give g\ = a + en while 
the third has the solutions: 

j^a + C12 + ci 3 a C7 , if c 7 / 0, 1 
g-2 = { c 6 a In a + C14 + ci 5 a , if c 7 = 1 
C6Q: + Cie In a + C17 , if C7 = . 

For the first possibility, we have finally the second case of 
Table 1 : 



9i = a 

32 = £a + , e # 0, 1 

33 = 1 + <5a + ^<5 a e 

where we have absorbed the constants cn,ci2 in the un- 
known function f , C11/4 + c 12 fs + fo — > /o, Eqs. (E(]), (|3|). 



After substituting these values of gi, <72,g 3 in Eqs. (|20| 
(P), we find that 



[/o-/6-/ 9 ]+ h+Us- h-Zh-&h a+ 



(46) 



/' 



/5/e- /s - (S0/9 



a E = 0, 



and 

P = {Po + Pia + P 2 a E ) 

= ff [/o + /i + (/4+C/5 + /2+C/3)« + M(^ + /3)a £ ] • 

By setting equal to zero the three expressions in the square 
brackets of Eq. (46) (since [i 7^ and 1 , a , a e are linearly 
independent vectors in the a-space for e 7^ 0, 1 ) we find 
the three R-relations for the functions G(R), M(R), fo(R) 
(which are the same with Eqs. (40), (45)). Using the func- 
tions fi,F and the definitions of Po and Pi we obtain five, 
first order , ordinary differential equations for G(R), F(R), 
M(R) and the two pressure components Pi(R) and Pq(R), 



dG 2 
dR 



dF _ F dAI _ F(F-2) 
dR 1-M' 2 dR 2R 



F 2 -4 2G 2 RP 1 
1-M 2 



(48) 



2fl(l-Af 2 ) 

P R Vi \(2M 2 - l) G 4 - M 4 + 2M 2 (l - G 2 )] 

M 2 (l-M 2 ) L\ J \ /I 



(49) 



dM 2 _ M 2 (l-M 2 ) ; 2 G 2 {l-M 2 ) 
— - (2M*-i)G*-M* \~ ed0V R 2 



^[{e + l)M 2 ~(e-l)G 4 }} 



+ 



(50) 



dR 



2R 2 G 2 + 2 ^G 2 (l 



1-G 2 



(l_ M 2) J 



dM 2 
dR 



M 2 F dF _ Su 2 _ M 2 (F 2 -4)(F-i) 
2R 2 G 2 dR R 2 M 2 4fl a G 2 ' 

(F-2)[(2M 2 -l)G 4 -M 4 ] 
^ i?G 2 M 2 (l-M 2 ) 2 



(51) 



dPg _ 2 dM 2 
dR ~ IS 1 dR 



v 2 _ 2M 2 (F-2) 
"WW 2 RG 1 



(52) 



Note that the third pressure component P 2 (R) is given ex- 
plicitly in terms of M and G (fs and /s). An integration 
of the above set of equations will give the complete solu- 
tion. However, this exercise is rather complicated since any 
physically accepted solution should pass through the various 
MHD critical points (Tsinganos et al 1996). This undertak- 
ing, together with a discussion of the solution and applica- 
tion to collimated outflows is the subject of the next paper. 
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Figure 1. Poloidal field and streamlines close to the stellar base 
for the asymptotically cylindrical 0-self similar model of case (2) 
from Table 1, for the following set of parameters: e = 0.5, v 2 = 
2SA4/r*V* 2 = 10, Si/ 2 = 3.5, 5 Q i/ 2 = 0.1, £ = -10, fi = 20, 
p* = (dM 2 /dR) = 1.6,F* = 1.1. 



It is worth mentioning at this point that our analy- 
sis of model (2) of Table 1 shows that mainly cylindrically 
collimated solutions are obtained. The set of Figures (1-2) 
illustrates such a typical solution for a representative set of 
the constants describing the particular model. This solution 
crosses the Alfven surface for appropriate values of the slope 
of the square of the Alfven number p* = (dM 2 / dti} , the 
expansion function P* and Pi* which satisfy the Alfven reg- 
ularity condition (Heyvaerts & Norman 1989, ST94) which 
is easily obtained from Eq. (A5) of Appendix A at (R = 
G = M = 1), i.e., 



P*P* = 2/ 4 



(53) 



The nonspherically symmetric part of the pressure Pi* is 
obtained from its definition while the functions /3* , /s* are 
calculated for R — 1 using the L'Hospital rule. Figs. (1,2,3) 
correspond to the set P* = 1.1 and p* = 1.6. Note that after 
the Alfven star-type critical point is crossed, the modified 
by self-similarity X-type fast critical point (Tsinganos et al 
1996) may be crossed by further adjusting appropriately the 
triplet of the variables (P* , p* , Pi*). It suffices to note that 



Figure 2. Poloidal field and streamlines as in Fig (1), but in an 
enlarged scale to show the asymptotical collimation reached after 
the oscillations have decayed. 



solutions crossing only the Alfven surface do not differ qual- 
itatively from those which in addition cross the modified by 
the present meridional selfsimilarity fast critical surface. 

Fig. (1) shows the shape of the streamlines on the 
poloidal plane and close to the Alfven surface. The cylindri- 
cal asymptotical shape of the poloidal streamlines is shown 
in the enlarged scale of Fig (2) . Note also the constant wave- 
length but the decaying with distance amplitude of the os- 
cillations, in full agreement with the analysis in Vlahakis 
& Tsinganos (1997). At the last shown fieldline ctout = 4, 
the toroidal fields vanish = 0, V$ = O.For a > a ou t , fi 2 
becomes negative , so there is no solution there. The same 
oscillatory behaviour can be seen in the fieldlines which are 
not rooted on the star but they are perpendicular to a thin 
disk around it (dotted curves in Figs 1, 2.) The oscillatory 
structure of all flow speeds before the flow reaches full cylin- 
drical collimation is also shown in Fig. (3) where we have 
plotted the characteristic velocities in units of the Alfven 
speed at the polar axis and Alfven sphere (a = ,R = 1), 
V*. 

The poloidal speed along the polar axis V PtPO i increases 
to a uniform superAlfvenic value and is higher than the same 
speed along the limiting streamline V Pt ii m (i.e., the last field- 
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40.0 



Figure 3. Outflow velocities in units of V*, the radial speed at 
the AlfVen point (a = 0, R = 1), for the parameters given in the 
caption of Fig. (1) of model (2) of Table 1. 



line rooted on the stellar base r star taken to be at 0.315r*). 
Both reach asymptotically uniform values after V Vt u m inter- 
sects the curve of the poloidal Alfven speed VA P ,um at R = 1. 
Note that corotation may be seen up to the Alfven distance 
R = 1: the azimuthal speed V^;j m at the 'limiting fieldine' 
increases until the Alfven surface is reached and drops from 
angular momentum conservation as the outflow expands al- 
most conically. Further away however, this speed too levels 
off to a constant value when full collimation is achieved, as 
expected. Finally, the fact that the jet has a large component 
of toroidal field is reflected by the large values of the Alfven 
speed associatd with the toroidal magnetic field, VA<p,Um, as 
compared to the rotational speed V^,n m . 



3 SYSTEMATIC CONSTRUCTION OF 

CLASSES OF RADIALLY SELFSIMILAR 
MHD OUTFLOWS 

To construct general classes of radially self-similar solutions, 
we make the following two key assumptions: 

(i) the Alfven Mach number M is solely a function of 0, 

M = M{0) , M(0*) = 1 , (54) 
and 

(ii) that the poloidal velocity and magnetic fields have a 
dipolar angular dependence, 

A=^A(a), a = -^- ) sin 2 0, R = ^ , (55) 

where B ,vj are constants. By choosing 0(9*) = 1 at the 
Alfven transition 0*, G(9) evidently measures the cylindrical 
distance vj to the polar axis of each fieldline labeled by a, 
normalized to its cylindrical distance vj a at the Alfven point, 



G (9) — vj/vj a . For a smooth crossing of the Alfven cone 
9 = 9+ [r = r a (a) ,0 = 0+}, the free integrals L and fl are 
related by 



^ = zul(A) = r 2 (a) sin 2 0* = zu 2 a . 



(56) 



Therefore, the second assumption is equivalent with the 
statement that at the Alfven conical surface, the cylindri- 
cal distance vj a of each magnetic flux surface a = const is 
simply proportional to \fa., exactly as in the previous merid- 
ionally self-similar case. 

Instead of using the three functions of a, [A,^a, fi) we 
found it more convenient to work with the three dimension- 
less functions of a, (qi , q2 , qi), 



da . 



qi (a)=J ±- 



92(Q) = 1| /" 2 *^a, 



q 3 (a) 



QM 



* 2 



l -da. 



(57) 



(58) 



(59) 



Following the same algorithm as in the previous case, we 
shall use (a , 0) as the independent variables and transform 
the derivatives with respect to r and to derivatives with 
respect to a and in the f- and ^-components of the mo- 
mentum equation. Integrating the resulting f-component of 
the momentum equation we get 

P = (h 1 aq 1 + h 2 aq 2 + h 3 q2 + h^qz + h 5 qi + h ) , (60) 



d2 , 
p = ££YPT 

8tt 



with 

P = [ ho h 5 hi h 3 hi hi ] , 
and 

Y = [ Yi Y 2 Y 3 Y 4 Y 5 Y 6 Y 7 } = 
1 qi aq 1 q 2 aq 2 q 3 aq 3 



(61) 



(62) 



and after substituting the pressure in the other component 
of the momentum equation we obtain 



Hh4.aq 3 + h 4 q 3 + h 3 (H - 2) aq 2 + h 3 q 2 + 



(^[i-Mf ) Qg ; 

(l-M 2 ) 



(63) 



+ h 5 qi + h = 0, 



where a prime in the functions of <Z;(a), i=l,2,3 and hi 

indicates a derivative with respect to their variables a 

and In sin 0, respectively, while the functions hj(0),j = 

1, 2, 3, 4, 5 and H are given in Appendix B. 

This expression is again of the form 

X 7 (0) Y T (a) + X 6 (0) Y 6 (a) + ■ ■ ■ + Xi (0) Yi (a) = , 



or , YX^ = 



(64) 
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with X the (1x7) matrix 
X = [ X\ X2 A3 X4 X5 Xq X-j 



hn h. 



(l-M2) 



h 3 h 3 (H - 2) h 4 Hhi 



(65) 



As in the previous case of meridionally selfsimilar so- 
lutions, we classify the various possibilities by the sets 
xxxxxxx. And, these sets always have "00" at the end, 
their first digit is "1", they have at most three "0's", while 
from the 2 7 possibilities we end up again with the cases 
1011100, 1101100, 1110100, 1111000, 1111100. Now the vec- 
tors qi(a), 92(a), 53(a) belong to a 3D a-space with basis 
vectors [ei(a) , 62(a) , 63(a)]. This space contains all vectors 
qi(a), i = 1,2,3 subject to the r-selfsimilarity constraint 
manifested by Eq. (63), i.e., that for a given such set qi(a), 
i — 1,2,3, the vectors l,aq i (a) , i = 1,2, 3 also belong to 
the same space. Each of the functions qi(a), i = 1, 2, 3 which 
satisfy this constraint are then a linear combination of the 
basis vectors ei(a), 62(a), 63(a). In the following, we choose 
ei = 1, e2 = <?i(a). All such sets of basis vectors give all 
possible radially selfsimilar solutions. Therefore, collecting 
all possibilities, we end up with the 6 classes of solutions 
shown in Table 3. 

In all of the cases of Table 3, from Eqs. @, (§|), (§i|) 
we find the form of the functions of a, 



.4 



2 QM q 2 - 
\l = — = — r a 



_ 2 B n zu 3 1 
*1 = ° a 2 q, . 
QM 



QMwJ^ah 
% 



(66) 



Finally, by substituting gi,g2,?3 in Eqs. ( p0[ ) , (63), we 
find the ordinary differential equations which the functions 
G{9) ,M(0) ,h (6) obey. 

In a-space, for each of the cases of Table 3 there exists a (3 
x 7) matrix K such that 

Y = [ ei e 2 e 3 ] K , (67) 
and from Eq. @ 

[ ei e 2 e 3 ] KX^ =0. 

If the basis vectors e, are linearly independent, then, 
KX^ = . 

These three equations are the ordinary differential equations 
for the functions of 9 in each model of Table 3, while for the 
pressure, 



P = fMd e 2 e 3 ]Kpt = ft (P 

87T o7T 



where 
Kpt = 



Piqi + Piez) 



[Pa Pi P2] 



t 



As with the previous meridionally selfsimilar solutions, the 
first two classes are of particular interest. The first corre- 
sponds to the following form of the free integrals: 



A = 



BoUj'i^/El £ , T ,2 CiBqZUo (F-3/2) 



o2 DtGM _ 
" = — 37^ a 



(68) 



This is a degenerate case, i.e., e3 = and we have 
only two conditions between the functions of 9. It fol- 
lows that we are free to impose a third relation between 
the unknown functions [G(0) , M(9) , ho(6)]. One possibil- 
ity is that such a third imposed relation is of the poly- 
tropic type, P oc p 7 (in this case ho = 0). In such a 
polytropic case which has been analysed in detail in Con- 
topoulos & Lovelace (1995), the magnetic flux is of the 
form A = f f (6)R F with f s {6) oc [sin 9/G(9)] F (for nota- 
tion see also Tsinganos et al 1996). The magnetic field at 



the equatorial plane 8 = 90° is B oc R 



the density 



p oc R 2F 3 , while the sound, Alfven and rotational speeds 
scale as their Keplerian counterparts, i.e., as _R _1//2 . Note 
that if [DiG(tt/2)/C*i] [(G 2 - M 2 )/G (l ~ M 2 )] 2 0=il = 1, 

the rotational velocity at the equatorial plane is exactly 
Keplerian. The classical and simplest subcase analysed in 
BP82 corresponds to the subclass with F — 3/4, wherein 
B oc 7?~ 5 / 4 . The two relations among the functions of 6 
are the two resulting first order differential equations for 
the Alfven number M{6) and dimensionless radius G(9) 
\m(x) = M 2 (6) and £(x) = G(0)/G(§) in the notation of 
BP82]. 

The second case is also degenerate since e$ — with again 
only two conditions between the functions of 6. As before, 
we are free to impose a third relation between the unknown 
functions [G(6) ,M(6) ,ho(0)], for example, a polytropic re- 
lationship. Then one can prove that this case is a subcase of 
the first one (if it is polytropic), for F — 2. All other cases 
shown in Table 3 are nondegenerate. 

The third class, is characterized first by a set of param- 
eters describing the particular model and the dependance 
of the free integrals on the magnetic flux function A (a), 
(xi , x% , Ei , Ei , G\ , G% , Di , Z?2). Second by the Alfven an- 
gle 9*. And third, by the set of the critical point parame- 
ters p* = (dM 2 /d8^J and ip+ which denote the slope of the 
Alfven number and the expansion angle, respectively, at the 
Alfven angle together with the pressure component Pi* 
through hs*. This triplet of 'dynamical' parameters fixes the 
physical solution and they are related through the Alfven 
regularity condition which is now obtained from Eq. (B9) of 
Appendix B at the Alfven angle 9* where M = G = 1 and 
^5 = i.e., . 



, tan(#* + V3*)p* 



(69) 



F 



QM 



As with the previous case of meridional selfsimilarity, this 
condition relates the slope of the square of the Alfven num- 
ber p t — (dA/ 2 /do) and the expansion angle 93* with the 
pressure component Pi* through Finally, the require- 
ment that the solution crosses the two slow and fast X-type 
critical points (modified by the radial self-similarity assump- 
tion, Tsinganos et al 1996) determines all these three 'dy- 
namical' parameters [<p+ ,p* , Pi*]. 

It is interesting to note that contrary to classes (l)-(2) in 
Table 3, this model (3) may be characterized by a scale, for 
example the radial distance on the plane of the disk where 
the magnitudes of the poloidal speed and magnetic field or 
the toroidal speed and magnetic field become zero. Hence, it 
occured to us that this is an interesting generalisation of the 
BP82 model and therefore worthy of further investigation. 
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Table 3. Radially Sclfsimilar Models 



Case qi(ce) 12(c) <73(«) constants 



(1) 


F-2" 


F-2" 


F-2" 


Ei , F - 2 ^ 




(2) 


Ei In a 


Di In « 


Ci lno 








(3) 


Eia xi + E 2 a X2 




Cia x i + C^a* 2 


e 2 + d'( - 


YC'i 


E 2 ,xi,x 2 ,xi — x 2 5^ 


(4) 


E\ In a + E 2 a x 


Di In a + D 2 a x 


Ci In a + C 2 cr I: 


E i + D 'i - 




x ^ 0,i = 1,2 


(5) 


Ei (In a) 2 + E 2 \na 


Di (In a) 2 + D 2 lna 


C*i (lna) 2 + C 2 lna 


e\ + d| - 






(6) 


Eia x ln« + E10F 


-Dia 31 In a + Z^a 1 


Cia 1 In a + C2Cr !: 


E\ + D\- 


YC\ 


7^0 



5.0 




equator 



Figure 4. Field and streamlines for the cylindrical r-self similar 
model of case (3) from Table 3 and the following set of parameters: 
xi = -0.6, x 2 = -0.5, Ei = -0.03, E 2 = 0.03, C*i = -1.5, C 2 = 
-0.6, Di = -25, D 2 = -10, 6»* = 62°, = 55°, p* = -3. At 
the disk level, V$ oc i? -1 / 2 while on the poloidal ficld/strcamlinc 
otout = 6.191736, Bp = V p = 0. 



Figs. (4-5) are a typical illustration of model (3) for de- 
scribing collimated jet-type outflows with an oscillatory be- 
haviour. In Fig. (4) the poloidal field and streamlines reach 
a cylindrical shape after undergoing oscillations in their ra- 
dius. As we move downstream, the amplitude of these os- 
cillations decays while their wavelength increases. In fact, 
the exact behaviour of the oscillations is analytically de- 
scribed in Vlahakis & Tsinganos (1997) where it is shown 
that they can be regarded as perturbations on an asymptot- 



40.0 



20.0 



0.0 



-10.0 



I I I 













v Ap 






I I I I 

0.1 1 10 100 



cotQ=z/m 



Figure 5. The characteristic velocities of model (3) of Table 
3 with cylindrical asymptotics are plotted in units of the z- 
component of the flow speed at the point (a = 1,8 = 7r/2) ,V 
and the same parameters as in Fig. (4). 



ically cylindrical shape which can be expressed in terms of 
the Legendre functions P^(cos8) and Q^(cos6>). According 
to this analysis, when fi 2 < 0, the asymptotically cylindrical 
shape is finally obtained through those oscillations. Then 
the perturbation (for f3 — > 0) is proportional to 6 )±M ~ I/ , or 
since a 2 < 0, proportional to " cos (| [i | In ^ + D )- 

In the example shown in Fig. (4-5) the amplitude of the os- 
cillations is rather weak. Note however, that cases also exist 
with an extremely strong oscillation amplitude and such ex- 
amples will be analysed in another connection. On the other 
hand, when y 2 ^ the asymptotically cylindrical shape is 
reached without such oscillations. Exactly this last possibil- 
ity is shown in the following case of Figs. (6-7). 

To further illustrate the various possibilities for the 
asymptotic behaviour of outflows starting from a Keplerian 
disk, we examine briefly the group of three models in Figs. 
(6-7), (8-9) and (10-11) where depending on the values of the 
model constants, we get one with cylindrical, parabolical, or 
conical terminal geometry: 

(1) In Figs. (6-7) a cylindrically collimated outflow 
(when 8 — > , (M 2 , G 2 ) — ► constants) is obtained for a set 
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Figure 6. Field and streamlines for the cylindrical r-self similar 
model of case (3) from Table 3 and the following set of parameters: 
X1 = -0.9, 22 = -0.6, Ei = -2.1421466, E 2 = 2.60994552, 
Ci = -3.2132198, C 2 = D 2 = 0, Di = -160.66099, 0* = 60°, 
= 74.704656°, = -1.1. At the disk level, V& <x R' 1 / 2 while 



on the poloidal field/streamline a ou t = 2, B p 



Vn = 0. 



of the model parameters: (xi, Ei, Ci, Di), i=l,2. The Alfven 
conical surface is taken at 9* — 60° where the slope of the 
square of the Alfven number is fixed as = —1.1 while the 
expansion angle y>* m 75° (the angle of the poloidal stream- 
line with the cylindrical radius). The characteristic scale of 
the model is taken to indicate approximately the radius of 
the jet, or more precisely, the distance along the disk where 
for ot ou t = 2 we have B p = V p = 0. In Fig. (7) the velocities 
on the reference line a = 1 are plotted in units of V , the 
z-component of the flow speed at the point (a = 1, 9 = 7r/2). 

(2) In Figs. (8- 9) an r-self similar model belonging to 
case (3) in Table 3 with parabolic asymptotical geometry 
(when 9 -> 0,{M 2 ,G 2 ) — > oo) is examined for another set 
of parameters (xi, Ei, Ci, Di), i=l,2. The Alfven conical 
surface is taken now at 9* — 45° where the slope of the 
square of the Alfven number is chosen as p* = —1.7 and the 
expansion angle 95* ~ 75°. 

(3) Finally, in Figs. (10-11) the r-self similar model of 
case (3) in Table 3 gives a conical asymptotical geometry 
for a third set of the parameters (an, Ei, d, Di), 1=1,2 and 




10000 



co<6=z/« 



Figure 7. The characteristic velocities of model (3) of Table 
3 with cylindrical asymptotics are plotted in units of the z- 
component of the flow speed at the point (a = 1,9 = tt/2) ,V 
and the same parameters as in Fig. (6). 



9* = 65°, = 75°, p 
exists only for 9 > 9 min where 9 m in 

r 2 



0.5. Note that now the solution 
17.5°. When this 

value of 9 is approached, {M 2, , G 2 ) — > oo. 

In all these four possibilities and along a given 
field/streamline, the outflow starts from the equator where 
V$ oc R -1 / 2 with a low subAlfvenic poloidal speed. This 
poloidal speed V p crosses the Alfven conical surface at 9+ 
in all cases. In the cylindrical case of Fig. (7), V p increases 
rapidly to a uniform value when collimation is achieved . The 
azimuthal speed V$ on the other hand, drops with height 
in all cases, as rotational energy is transformed to poloidal 
kinetic energy. Finally, the azimuthal Alfven speed is the 
strongest in the cylindrical case where the toroidal magnetic 
field is responsible for the ensuing final collimation. 



4 SUMMARY 

In this paper we have examined a systematic way for con- 
structing exact MHD solutions for plasma flows. The first as- 
sumption was to consider the ideal plasmas MHD equations 
for time-independent conditions, Eq. (^-^), without impos- 
ing the extra constraint of the frequently used polytropic 
assumption. Second, we confined our attention to axisym- 
metric situations in which case the poloidal magnetic and 
velocity fields can be expressed in terms of the magnetic 
flux function A while several integrals exist, Eq. (^-^)- In 
that case, besides A, a second natural variable is the Alfven 
Mach number M, Eq. (^). We denoted by G the cylindri- 
cal distance w of a poloidal streamline from the system's 
symmetry axis, in units of the cylindrical distance of the 
Alfven surface from the same axis, zu a . Third, we further 
confined our attention to transAlfvenic outflows in which 
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equator 



Figure 8. Poloidal field and streamlines for the parabolic r-self 
similar model of case (3), Table 3 and the following set of param- 
eters: xi = -0.9, x 2 = -0.6, Ei = -0.8252542, E 2 = 1.00547, 
Ci = -1.23788, C 2 = D 2 = 0, Di = -12.378813, 6>* = 45°, 
</3* = 75.465545°, p* = -1.7. In this case Vj, <x i?" 1 / 2 on the 
equatorial plane while on the streamline a ut = 2, B p = V p = 0. 



case the regularization of the azimuthal components in Eq. 
(§~H) requires that the ratio of the two integrals of the to- 
tal specific angular momentum in the flow L(A) and coro- 
tation frequency Q(A) is some function a(A) [as in Eq. 
(@)]. By introducing some reference scale zu this function 
q is dimensionless, [as in Eq. (Q) where zu = r*]. Ap- 
parently (M, a) is a rather convenient set of dimension- 
less variables for describing all physical quantities in the 
poloidal plane. For any set of orthogonal curvilinear coordi- 
nates suitable for describing axisymmetric problems, we may 
then convert their poloidal coordinates to (M,a). Exam- 
ples are, spherical coordinates [r(M, a), 6(M, a),<j>], cylindri- 
cal coordinates [z(M, a), w(M, a), 0], toroidal coordinates 
[u(M,a),v(M,a),4>], oblate/prolate spheroidal coordinates 
[£(M, a), T)(M, a), <f>], paraboloidal coordinates, etc. Then, 
the distance from the symmetry axis of the outflow is 
G(M, a). In the present first study we made the simplifying 
fourth assumption that G is independent of a, G = G(M) 
only. Finally, to re-establish the connection with the geom- 
etry of the problem and the particular set of the coordi- 
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Figure 9. The characteristic velocities of model (3) of Table 
3 with paraboloidal asymptotics are plotted in units of the z- 
component of the flow speed at the point (a = 1,0 = 7r/2) , V 
and the same parameters as in Fig. (8). 



nates used, we made our fifth and final assumption that 
M = M(x) (and G = G(x)), where x = r , or , X = 
This leads then to the two broad classes of meridionally and 
radially self-similar outflows. Needless to say that additional 
symmetries may in principle be considered, something which 
may be taken up in another connection (equilibria in toka- 
mak geometries, etc). 

After these five assumptions are well posed and with 
the help of a simple theorem, it is possible to (i) unify all 
existing exact solutions for astrophysical outflows (Tables 
1,2 and 3) and (ii), to qualitatively sketch a few of them. 
With this method, the system of the coupled MHD equations 
reduces to a set of five ordinary differential equations for 
the dimensionless jet radius (G), the flow's expansion factor 
or angle (F, or ip), the Alfven Mach number (M) and the 
two pressure components (Pi and Po)- The requirement that 
the solutions pass through the Alfven critical point gives 
a condition relating the values of the expansion function 
or angle, Alfven number slope and pressure component at 
this critical point. The Alfven regularity conditions, Eqs. 
( p9[ ) , ( j53| ) is similar to that discussed in Heyvarts & Norman 
(1989) and ST94. 

As a byproduct of this construction, two representative 
models for radially and meridionally self-similar outflows, 
BP82 and ST94, respectively, have been generalized. In the 
former case of BP82, it is well known that the cold plasma 
solution is terminated at a finite height above the disk while 
the general case (3) in Table 3 extends all the way to infinity. 
Also, it is shown that the expressions of the MHD integrals 
which correspond to the ST94 model are only a special case 
of case (2) in Table 1. 

Having in mind the ubiquitously observed collimated 
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Figure 11. The characteristic velocities of model (3) of Table 
3 with conical asymptotical geometry are plotted in units of the 
z-component of the flow speed at the point (a = 1, 9 = tt/2) , V 
and the same parameters as in Fig. (10). 



sure, radiative forces, etc). For example, the same formalism 
may be used also in radiation driven winds. 



Figure 10. Field and streamlines for the conical r-self similar 
model of case (3), Table 3 and the following set of parameters: 
xi = -0.1, x 2 = 0.01 and Ei = -78.601635, E 2 = -728.31337, 
Ci = -4.3231, C 2 = D 2 = 0, Di = -43.231, 6»* = 65°, 
ip* = 75.784234°, p* = -0.5. In this case oc R' 1 / 2 on the 
equatorial plane while on the poloidal field/streamline a ou t = 2, 
Bp = V p = 0. For large distances from the disk all lines with 
a > go asymptotically to the line 9 = 9 m i n . 



outflows from astrophysical objects, we paid more atten- 
tion to the selfconsistently derived asymptotical shape of 
the streamlines. Of the various such asymptotic geometries 
derived, a prominent member seem to be the cylindrically 
collimated jet-type solutions, in accordance also with the 
conclusions of observations (Livio 1997) , general theoretical 
arguments (Heyvaerts & Norman 1989) and recent numeri- 
cal simulations (Goodson et al 1997). Another feature that 
appeared in the solutions is that cylindrical collimation may 
or may not be achieved with oscillations in the width of the 
jet (Vlahakis & Tsinganos 1997). Although in the examples 
analyzed here the amplitude of the oscillations is rather weak 
and the flow collimates rather smoothly, preliminary results 
show that cases also exist where it can become rather large 
and the final radius of the jet can be much smaller than the 
initial large cylindrical radius and corresponding opening 
angle. Finally, we should note that the pressure P denotes 
the total pressure (including gas pressure, Alfven waves pres- 
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APPENDIX A: FUNCTIONS OF R 
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sin 6 sin (tp + 6) 
cos 8 cos (ip + 9) 
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cos 6 cos (ip + 9) 
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where the expansion angle ip is the angle between the line 
and the equatorial plane, which is a function of 9. 
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